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ON  SHOCKS  ATTACHED  TO  A  WALL  IN  A  STATIONARY 
OR  PSEUDO -STATIONARY  TWO-DIMENSIONAL  FLOW 

by 

Y„  So  Chow 


L.   Introduction 

Shock  waves  attached  to  a  solid  "boundary  have  "been  discussed  by  many 
authors  in  stationary  or  pseudo-stationary  two-dimensional  flows,  In  the  stationary 
^ase,  Tsien  [9]  showed  that  at  the  attachment  pointy  0,  of  a  normal  shock,,  the 
Kach  numbers  before  and  after  the  shock  are  determined  by  the  ratio  of  the  curva- 
tures of  the  boundary o  He  assumed  that  the  flow  ahead  of  the  shock  is  isentropic 
and  irrotational  and  that  all  the  quantities  involved  are  finite  and  continuous 
except  in  crossing  the  shock,  Lin  and  Rubinov  [3]  extended  his  result  by  allowing 
the  curvature  of  the  shock  and  the  pressure  gradient  behind  the  shock  to  become 
infinite  at  the  point  0.  Let  M  =  ~~  +  —  y(y+l)^  +  k   ,  They  proved  that  over 
a  convex  (concave)  boundary  with  continuous  curvature  a  normal  shock  can  be  formed 
only  when  the  Mach  number  ahead  of  the  shock  is  above  (below)  M  „  By  assuming 
finiteness  of  the  curvatures  of  the  boundary  before  and  after  the  shocks  they 
found  that  the  normal  shock  in  a  neighborhood  of  the  point  0  is  approximately 
x  =  Py  '  ,  where  p  is  a  constant,  y  =  0  is  the  tangent  of  the  boundary,  x=  0  is  that  of 
the  shock,  and  0  is  the  origin,,  At  a  slightly  earlier  time  than  Tsien 's  paper, 
Thomas  [5]  obtained  a  relation  between  the  curvatures  of  the  stream  lines  behind 
the  shock  and  those  of  the  shock,  and  later  he  [6]  extended  the  result  to  the 

derivatives  of  the  curvatures ,  if  the  flow  is  stationary  and  the  flow  ahead  of 

u 

the  shock  is  uniform,  Taub  [k]   gave  a  simpler  derivation  of  Thomas'  results  and 
extended  them  to  pseudo-stationary  flows. 

Brown  [lj,  using  devices  due  to  Thomas,  extended  Lin  and  Rubinov* s 
results  to  stationary,  iso-energetic  flows.  The  author  extends  Taub's  treatment 
to  obtain  a  general  relation  among  the  curvatures  of  the  streak  lines  behind  the 
shock,  before  the  shock,  and  of  the  shock,  for  a  general  stationary  or  pseudo- 
stationary  flow. 

Let  the  boundary  be  a  straight  wall  and  let  the  flow  ahead  of  the  shock 
be  uniform.   In  part  k,   the  author  proves  that,  if  the  shock  is  normal  to  the  wall 


r 


at  0,  the  derivatives  of  odd  order  of  the  curvature  of  the  shock  are  zero  at  0, 
provided  that  the  shock  is  regular  there.   In  part  5,  the  curvature  of  the  shock  is 
allowed  to  he come  infinite  at  0  and  approximations  for  the  streak  lines  and  the 
shock  near  the  point  0  are  obtained. 

2.  Main  Equations 

(A)     Let  Juhe  a  shock  in  a  two-dimensional  stationary  or  pseudo-stationary 
flow  of  an  ideal  gas  with  constant   specific  heats  and  without  heat   conductivity 
and  viscosity,    (x..,xp)   a  fixed  cartesian  coordinate   system  in  the  plane,,    and  the 
shock  he  described  by  the  parametric  equations 

x±     =      (1  -  e)a.(s)  +  st  ai(s)  (2.1) 

where  i  =  1.2,  s  is  the  arc  length  along  the  shock,  and  e  =  1  or  0  according  to 
whether  the  flow  is  pseudo-stationary  or  stationary.  For  simplicity,,  the  summation 
convention  will  he  used,  i.e0,  in  a  term  repeated  indices  are  dummy  indices  of  the 
summation  from  1  to  2,  and  all  the  free  indices  may  take  values  1  and  2  except 
when  the  contrary  is  mentioned.  We  will  assume  that  continuity  and  differentiability 
requirements  on  the  functions  involved  in  the  discussion  of  this  section  are  always 
satisfied. 

Choose  independent  variables  (s,T)  by 

•-■'  x±(s,T)  «  (1  -  s)a1(s?T)  +  €tai(s,T)  ,  '(2.2) 


where 


aa.(s,Tj     U 

X 


for  I  ^  0,  U^SjT)  -  mi(s,T)  -  ea±(s,T),  (2.3) 

a,  (s,0)  =  a,  (s)  (2.*0 


1 


and  u.(s.,T)  and  u.  (s,  =  T),  for  T  >  0,  are  the  velocity  components  relative  to  the 
fixed  coordinate  system  (x_  ,x  )  after  the  shock  and  before  the  shock,  respectively. 


Define 


3a.(s,T)              \ 

Ms,T)  =  ■ — ^ -,  *-(s,T)  =  /\,  K    3   tan  a(s,T)  =  r~   (2.5) 

U,(s,T)  =  X.  U.   and  U  (s,T)  =   V  U.  ,                  (2.6) 

tv  9    '             ix  nx  '  11* 


_2= 


■where 

-P    (s,T)  =  \     ,     n)    (s,T)  =  -X.    .  (2.7) 


1 
Put 


tan  |i(s,T)  =  —   and  tan  J2(s,T)  =  ™  .  (2.8) 


Then 


^  2(s,T)  =  a  +  u  -  -—  .  (2.9) 

In  the  following,  partial  differentiation  with  respect  to  s  and  T 
will  he  denoted  "by  a  subscript  after  a  comma,  for  example,, 

3, 


t        m      SS.  f  2  -         2    f 

•S      9S        *S  X     a2s3T 


The  curves  T  =  constant  in  the  (a,  ,  a  )  plane  will  he  called  displaced 
shocks  and  s  =  constant  streak  lines „   It  is  obvious  that  the  shock  itself  is 
the  curve,  T  =  0,  and  in  the  stationary  case,  streak  lines  are  the  stream  lines. 
We  will  put 

k(s,T)  =  c^s   and   K(s,T)  =   ^  2,T  or  "^2,1  (2ol0) 

according  to  whether  U  (s,T)  >  0  or  IT  (s^T)  <  0„ 

(B)  Fundamental  equations.  The  hydrodynamical  equations  for  the  flow- 


are 


dip**) 

TT-   +2£P   "    Q    '  (2°ll} 


i 


/>(U.   ~~+  €  U.)   +4f-     ■     °     »  (2.12) 

Ui  4i  -   S,T   -   °   '  (2-i3) 

where  /^p  and  S  are  the  density,  pressure  and  specific  entropy,  respectively.  The 
Sankine-Hugoniot  equations  (along  the  shock)  are 


1 
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Ut(s,0)   =  U£   , 

tun(s,o)=u-  , 


and 


y(s,o)  -Jr-1  ^(cr2.!)  , 


7^(8,0).  4=  -^^ 


P         2+(7-l) 


(T 


(2.1^) 
(2.15) 

(2.16) 
(2.17) 


where  f(s,0)  and  f'(s,0)  of  the  variable  f  are  defined  by 

f(s,0)  =  11m  f(s,T)  ,     ff(s,0)  =  lim  f(s,T)  , 

U* 
7  is  the  ratio  of  specific  heats,  o-(s,0)  =  -2-  =  mf  cos  u',  m(s,T)  ■  -  and 

c2(s,t)  -  a. . 


Let 


Ji^T)  =  i  log  p,   Cl(s,T)  .  0,   C2(s,T)  =  ^22L 


and 


From  (2.11)  to  (2.17),  ve  have  [4; (5.10)] 


(2.18) 


where 


\(l-m  cos  |i; 


sin  u 


m  cos  p, 


(l  -  -^-)cos  u   sin  u. 
m 


(2.19) 


IK^T)!)  -  [|J^(.,T)||^-X 


sin  (i.  -m  cos  (j. 

(-jp  -  l)  cos  n   sin  [i 
m 


(c)  Derivation  of  the  main  equations.  From  (2.l8),  along  the  shock, 

*],    m(s,0)  =  A.  .  5.  (2.20) 


^i,T(3,o)=A:.  5' 


(2.21) 


-U. 


By  differentiating  (2.l6), 


h(T(r- 


51(s,0)    -  a{(s,0)   = 

27  cr     +1-7 


=  B1,    say. 


(2.22) 


From   (2.1*0   and   (2.15), 

tan  |i(s,0)   =  -7)  tan  \i} 

and  then 

52   sec  n(s,0)    -  ^6^sec^'(s,0)  =  k(sec^   ~  77  sec'V  )   +   7)        tan  n' 

(2.23) 


A.i 


?  /        2                     2    ,\            ^  tan  u  cf  s 

=  k(sec  |i   -  T)  sec  (i'J   +  — =L_  „   — ^—5 

2+(7-l)<r2 


a— 


S2,  say, 


(2.20)  to  (2.23)  give  six  equations  for  eight  variables  \.      9   1?    5  and  5.'. 

^s  XoX        ^  X  «  X  1  X 

Therefore  ve   can  solve   for  5.    in  terms  of  ^  !      ,   k,   B  _,   B     and  the   flow  variables, 
The  matrix  of  coefficients  of   (2.20)   to   (2.23)    is 


Jl,T  jfg,T 


5i 


52 


-All 

-A12 

1 

"A21 

-A 

A22 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

2 
sec  (j. 

0 

0 

0 

0' 

1 

0 

0 

0 

A' 

All 

A' 
A12 

0 

A1 

A21 

A! 
A22 

0 

-1 

0 

0 

0 

-1) 

0 
0 


1,T 


2, 


^sec  \i 


a, 


B. 


Let  D  be  the  determinant  consisting  of  the  first  six  columns  of  this 
matrix,  and  D  be  that  obtained  by  replacing  the  first  column  in  D  by  the  last 
column  of  this  matrix.  Then 


D(s,0)   = 


sec  [i 


■r  1  2  ■        2    , 

1   -  m'      cos  (j. ' 


Dl(s,o) .  (a^2  jijT  -  Aia  y2. 


B, 


2         2 

1   -  m'      cos  p.! 


N  2 

-jsec  (j. 


• 
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Therefore, 


51(s,0)  =  sin  u'    .    lj^T   -  mV  J^T  +  B±   .  (2.24) 

Similarly,  2  2 

B2(s,0)   «^£2|it(-i^ilCOBll..   J.      )   +Bpcos%    .  (2.25) 


cos  |i  ■       m' 


Let 

0(8,T)   =  U1)2    -  U2>1        . 

Then,    [4; (4.19)] 

Un0(s,T)  =  W  s   +  -^  +  e  Ut      .  (2.26) 


Let 


Then 


2 
h(s,T)  =  \-  +  —5-  c2     . 

o  v         v  p 

h    (s,t)  =vv     +^S-C      =  rr     +  -M-*^  -  —4)  (2.27) 

,sv    '    '  9s       7-1      ,  s  ;s        7=1     />  /p  2'  v^.t.i  t 


s        7=1    ^  ^ 


Therefore 


y  c  v  ,     ,  w2 

^  v  c    ){S>1)  ~     v   L1  +     0  2        J    "       2     h,s 

2c  2c  ' 


Hence 


where 


.  K (Un  t-  If  -  e  ut)[i  +  iZi^d]   -  Zr|  h      . 

v        n  /^  t  2c2  2c2     >S 


v  c  €U 

(-£  -  -^)(s,T)  =    Aj  +  f±)Vj   -    (ig  +  Z*>  j  -  -g*     B  (2.28) 

v  2c  m  J  v 


B3(s,T)  =  2^%(h  h  +  e  U+)    .  (2.29) 

Now  [2;   p. 10] 


,3v,,,,   -  2c2^s   T  t   ut 


where  c     is  the   specific  heat   at  constant  volume.      Hence, 

S      (S,T)  =   c    (^2.  .  1    f?    )    .  (2.30) 

,s''?/vNp^/,s/ 


.6- 


From   (2.26)   and   (2.27), 

Un0(s,T)   =    fc.   +  £Ut)    -  T7-E1_  ( 


Hence, 


=  7^T(B3  -a7T-s,s)- 


B,(8,T)-£|V  +     1     s>s 


3 


(2.31) 


Since  cr  (s,0)  =  m1    cos  li',    we  have 


«r-f8(B,0)  =  m;s   cos  ^    -  mVs   sin  tf  =   (-*»  -  _^)[r  -(j'^-kV    sin  ^    . 
By  (2.28)  and   (2.31), 


1  7-1 


crB(s,0)  =   -(—     +  2^-)cr  5|   -  m'    sin  |a '   5^  +  km'    sin  \x 


/COS    LI  '       /,  1  , 

+  (    „,      0'  -  ^r —  s1  )o-  . 


V 


27c  ,S' 

V 


Nov  [1*;(5.7)] 


where 


B.(s,T)   =  a.  .   \.  _ 
iv    '  ij  Jj,T 


||a      (s,T)||    =  \[ I 


sin  li 


-m     cos  li 


(-^  -  l)cos  [i        sin  li 


m 


(2.32) 


(2.33) 


Thus 


<T,(»,0)   =   (^-  -^)      sin  n'    j^T  +   (cos  2^   +  2^1  ^2^,  j,^ 


(2.3*0 


+  ta'    sin  ll'   +  <r(£2*J^0  .  -±-_  s,    } 


27c        ,  s 
v 


From  (2.2^),    (2.25),    (2.22),    (2.23)   and   (2.34),   ve  have 

5,(s,0)  =  d..\\  _  +  e.k+  f.CSSSJil^.    .      1      S.    ) 
1  iJ   oj,T  1  iv     v'        r  27c        ,s'' 


(2.35) 


where 


-7- 


U 


sin  \i 


,   6d-  +  8  cos  [i'+  I-7 


27  or  +  1-7 


2                   2 
ij-  cos  u.1    -  2CT     -I+7      , 
c *-  m'o— 

27  ~    -  1  -  X 


sin  2|i       C    1 


2   sin  n1 , 

'    Til 


-    1    + 


[(2-y)m,2-k]\,  sin2^^cos2^'^3^1>r-2+2] 


[2+(<r-l)a-2]m'2 


2  cos  pi'    [2+(ar-l)a-   1 
(2.36) 


and 


el       fl 
e2       f2 


^m1    sin  |i' 


27  o~  +1-7 


1  - 


T( 


2 

COS    [I 


2  sin     \i 


2  2 

cos  [i1        2+(iT-l)a- 


m'    sin  uj 

cr- 


ko-2 


2j(f    +1-7 


2  sin  2\i 

2+(r-i)a_; 


2        2 
■ — —  m     sin  2p. 
7+1 


i^o-2 


7+1 


2*^    +i-r 

2  sin  2[i 
2+fr-l)^2 


(2.37) 


Differentiating    (2.35),   ve  have 

5^8n(s,0)  =  ^^,snT 


■where 


(2.38) 


sJ'W)  =  ^  +  4ia,S  Si.^  +  %s*,s»-i+<^-) f8  0;^ 


and 


f . 

i*s  S' 
27c        ,sn 

7    v     ' 


;[0)(s,o)  =  0  . 


(2.39) 


(2.1*0) 
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Now   [4;    (7.5)   and   (7.6)] 


(2.41) 


where 


P.      vs,Tj  -  P^       +Aij^  5.^n.2+A..  ^    (   r    )A.k^r  B^n-l- 


■eA^"1)   C.  ^.2 


and 


Pr(s,T)   =  0 


(2.42) 


(2.43) 


Since 


]j,sVS>°>  -  5j,s-lT  +  £Cj,Sn-lT  "  -S'tfi,*-*   -Pin))+^,«n-lT)    (2.W) 


5  (s  0)  -  d     a'^T'  i    +  p  k         +  f   r^os  ^V  1  .     q,        -  Uh^n^ 


where 


h.      <s,Oj  -  g„  d^a^  P^       -€C^sn.lT;    . 


(2.46) 


From  (2.45)   and   (2.4l),   we  have  the  main  equations 
I.   __(s,0)   =  A?n)d.    a'(^l}  >;   mn  +  A?n)e.k     n.x   +  Afn)f.(^_M^.  .  _J_  rf.      ) 

(2.47) 


(n 


where 


(n)/      ^\         «(n)    1  (n_l)        t,(q)  /^   I  ON 

<y    y(s,0)   =  A^i  '   h^        '   +  P£    '    .  (2.48) 

A  polynomial  r(s,0)  of  \,^l,   k  B£-l/  0  s£-l  and  S  g£,  for /l=  1,  2,  ...,  n 

without  constant  term  and  whose  coefficients  are  functions  of  flow  variables  will 
be  denoted  by 

r(s,0)  =  r(  (   jl  ,  k  o_i,  $     fl.i,  S  0   ,   1  <  jl<  n)  .       (2. 49) 
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Then 

b^coj  =  41)(^,s> k>  ^  s;s)  .  (2.50) 

since 

and  by   (2.28) 

Since    [It;    (5.12)  to   (5.15)] 

I  '*,*<■'*> -TTcji,s  (2-5l) 

J  M,T(S'T)  "   "(i  +  TT  m)Jl,T   -  I     «  (2-52) 

\      (s,T)  =  cos  (lj2  a      ,  (2.53) 

hj1)    (s,0)  =  41}    rfj^,   k,    0',    B«a)  (2.55) 

Pp}(s,0)   =p(2)(j.^    k,    0-,    S6)    .  (2.56) 

Therefore 

<42)(s,0)  =  <£2)(fctS,  k,   0',   S's)    .  (2.57) 

An  induction  argument   shows  that 

q^n)(s,0)   =   q[n)    (fajL-1,    k,s^2>    0y-2,    ^V"1'    2</<n)    .      (2.58) 

Equations  (2.V7)  extend  the  previous  results  obtained  by  Thomas  [6],  Brown  [l^ 
and  Taub  [h]   to  the  general  stationary  or  pseudo-stationary  flows.  The  coefficients 


and 


we  have 


and 
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(n)  (n) 

AK    J   e  have  been  denoted  as  A^   in  [k;    (10.9)  and  (10.10)]  and  A^e.  as  G   (m»  fa) 

"  £-j   i     n-1 

in  [6;  (19)],  provided  that  the  flow  ahead  of  the  shock  is  uniform,  i.e., 

Ul(s,T)  =  constant  and  u2(s,T)  =  constant  for  T  <  0.   In  [7],  Thomas  showed  that 
for  m(s,0)  >  1 

A^V^O)  =  Gn-1(m,,a)  5/  0       for  n  ^  1,2,  ...,  n,      (2.59) 

but  if  a  =  a±   is  such  that  m(s,0)  <  1,  then  any  neighborhood  of  a  contains  a 
value  of  a  at  which 

Gn_1(mSa)  *  0  (2.60) 

for  some  positive  integer  n  . 

(D)  Simple  consequences  of  (2.47).   Let  the  flow  ahead  of  the  shock  be 
uniform.   Then 

ji,T(s,T)  =  0  (2.6l) 

and  (2.47)  for  i  =  2  becomes 

j2,Tn(s,0)  =  A^e.  k^n=1  +  4n)  .  (2.62) 

If  m(0,0)  >  1  and  K(0,T)  =  0  for  T  >  0;  then  by  (2.59);  (2.62)  and  (2.58) 

k(0,0)  =  q^2)  .  k^  =  q^3)  =  k^2  =  ...  =  k^n  »  ...  =  0  .    (2.63) 

It  is  easy  to  see  that  (2.63)  still  holds  if  we  assume  only  that  the  continuity  and 
differentiability  requirements  on  the  functions  involved  in  obtaining  (2.^7)  for 
s  >  0  are  satisfied  and 

lim  m(s,0)  >  1,   lim  K  ^(8,0)  =  lira  K'   (s,0)  »  0  , 
s^O  >  }l 

limJl  Tn+l(s>0)  =  li™  0,sn(s>°)  =  llm  S'sn+l(s,0)  =  0  and 

lim  k  n(s;0)  =  k  n(0,0)    for  n  =  0,  1,  ...,  n 
?  s  j  s 

If  the  flow  ahead  of  the  shock  is  uniform  and  after  the  shock  supersonic  near  the 
vertex,  then  (2.63)  gives  proofs  to  the  following  results  at  the  vertex,  (a)  The 
attached  shock  to  a  cone  at  its  vertex  is  itself  a  cone,   (b)  The  attached  shock 
to  a  wedge  is  straight.   (c)  In  regular  reflection,  the  reflected  shock  behind  a 
straight  incident  shock  is  straight. 
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From   (2.17), 

p  (s.o)  = : i±—— W/°sjd  ii'    , 

By   (2.30)   and   (2.1*5), 

S,s(s>°)  -  »T(7  JljS  -  ^)  -  V^e,  -    .  »»•  f  °  |_)   „         (2.65) 

'  (2+70--  o-^)c^ 

If,   moreover,   0(s,O)  =  0,    then  by   (2.1*5) 

f ' 
0  =  6!(s,0)  -  d '     a       T         +  e»k  i     Q 


v 
Therefore 


f±  2f„'m'    sin  u< 

°  =    (dljej+ei   ~  T  el  +  T~T~ 2— )^^°)   =  "  ,    "ay. 

(2+/T-  -  0-^)(T"  X 


(2.66) 


(2.67) 


X.  are  functions  of  m'(s,0)  and  sin  h'(b,0),  whose  coefficients  are  constants. 
lince  the  Jacobian  of  X±   and  Xg,  with  respect  to  m«  and  sin  u',  is  not  zero  as 
1  -  0,  Xx  and  X2  do  not  have  a  relation  between  them.   For  each  interval  I,  if 
X^O)  ^  0  or  X2MV0forsomes0cI,  then  k(S(),0)  =  0  by  (2.67),  and  if 
X1(s,0)  =  X2  =  0  for  all  s  6  I,  then  m'(s,0)  .  constant  and  u'(s,0)  =   constant 
'  s  e  I  and  the  flow  is  stationary  and  k(s,0)  =  0  for  s  €  I.  Therefore,  for 
ch  interval  I,  there  always  exists  an  sQ  e  I  such  that  k(S(,0)=  0.  By  the 
continuity  of  k,  we  prove  that 

k(s,0)  =  0        for  all  s  ,  (2.68) 

Thus,  if  the  flow  ahead  of  the  shock  is  uniform  and  the  flow  after  shock 
irrotational,  then  the  shock  is  straight  and  the  flow  after  the  shock  is  uniform 
also.  This  result  can  also  be  obtained  from  a  much  stronger  result  of  Taub  [Univer- 

■ty  of  Illinois  lecture  notes].  Every  irrotational  non-isentropic  pseudo-stationary 
flow  has  constant  particle  velocity.   But  the  proof  given  here  is  simpler. 

3.  Normal  Shocks 

In  this  section  we  are  going  to  discuss  a  normal  shock  attached  to  a 
)lid  boundary  with  continuous  tangent.  The  curve  s  =  0  will  be  the  boundary 
diich  has  continuous  curvature  everywhere  except  at  (0,0)  and  there  it  may  have 
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a  finite  jump.  The  shock  a.(s,0),  s  >  0,  is  regular  everywhere  except,  possibly, 

at  (0,0),  where  lim  k(s,0)  =  k(0,0)  may  be  infinite.  The  flow  in  front  of  the 
s— >0 

shock  is  regular  off  the  boundary.  We  assume  that 

lim  f(s,0)  =  f(0,0)   and   lim  f ' (s,0)  =  f • (0,0)  (3.1) 

for  all  quantities  f  used  below,  and  that 

0'(O,O),  S«s(0,0),  K'(0,0)  andp»T(0,0)  are  finite.  (3.2) 

Since  the  shock  is  normal  and  the  boundary  is  continuous  at  (0,0), 

a(0,0)  =  %-,    sinu(0,0)  =  sin  ;i.J(0,0)  =  0„  (3.3) 

(A)  Limiting  case  of  (S.lff)    for  n  *  1.     When  n  =  1  and  i  =  2,    (2.^7)  be- 
comes 


lim  Ad     (s,0)  =  0, 


6*0 

li.  A2J  i    (.,o>  -  gVg  -3-r?   (0;0, 

m   (27m"   +1-7) 


2.,n  „,„  2 


lim  Ag     e  k(s,0)  =  .-  &!U£  llffi  k  sin  ^  _2C^j^  iim(^       ,-      )sin  ,, 
^0     ^      J  (7+i)m2  (7+l)m2v  J 


2v'f2m2+l)   ,4    f  .,  v,2f2m2+lj  ,  ,    ,       2  ,v 

=  -  ^-    ^  j_  lunfcos  pi9)        «  p      ^  limfcos  u'/ 

(7+1  )m  v  ^s      (r+J> vF'n  * 


2 
lim  A_„f  .(s.O)   =   -  ifm'    Y'        _/q  o'. 
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From   (3.4),    (2.8),    (2.9)   and   (2.10)   we  have 

K(0,0)   =  _il [(2m.2-3-7)K.    -Ag.        .  I&J    (0  0) 

nT^/m'   +I-7)  7Cy       ,a         yI2     J    ^0'0> 

v'(2m2+l)    ,  p 

-/  o        (0>°)  lim(cosdn')      (s,0) 


{nlJA      ^"^"^.^    .  (3.5) 

(3.6) 


(—  +  7   -  1)K(0,0)  >      [(2m'2   -  3   -  7)K-    --2-S'        -J_*](00) 
m'  re  ,s  2     J  vw->u' 


(3.5)   and   (3.6)   hold  for  a  general  stationary  or  pseudo -stationary  flow. 
The  equality  of   (3.6)   holds   if  and  only  if 

lim  k  sin  u    (s,0)   =  0      .  ,_  _v 

s^O  ^3.7; 

(B)  Brown's  result.  Let  the  flow  he  stationary  and  iso-energetic.  Then 
by  (2.29)  and  (2. 31) 

7rS',s  =  -H"'J'   >  (3.8) 

and   (3.6)  reduces  to 

(—  +  7   -  1)K(0,0)  >    (2m'2   -  3   -  7)K'(0,0)    -  2(*£  +  -%0'U'    (0,0) 
m  c2         v'2  n 

(3.9) 
Put 


r  :  v,"tf     ■  (0,0)   . 
2,T 


T'Y 


Then   (3.9)  becomes 
2 


(~2  +  7   -  1)K(0,0)  >   (2m«2   -  3   -  r)K'   +  2[(7-l)m'2  +  2]/V    .      (3.10) 

m'  v  ' 

Equation    (3. 10)   is  the  same  as  Brown's   [l|(a6)J    ,   which  should  he 

2 

K(0,0)  >~   m  (M2   -  ^±3)K'(0,0)   +  2(y-l)  M2  K«      (0,0)    .         (3.ll) 

(*-l)M+2  d 
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(C)  Suppose  that  the  flow  is  irrotational  isentropic  ahead  of  the  shock. 
Then  (3.6)  reduces  to 


(—  +  7   -   1)K(0,0)  >  (2m'2  -  3  -  7)K'(0,0),  (3.12) 


which  extends  Lin  and  Rubinov's  result  ["3;  (^.7)]  to  pseudo-stationary  flows.   If, 
moreover,  k(0,0)  is  finite,  then  (3.7)  is  satisfied  and 


2 


(—  +  7  -  1)K(0,0)  =  (2m' d  -  3  -  7)K'  (0,0),  (3. 13) 

m' 

which  was  proved  by  Tsien  [9',(l6)]   and  Lin  and  Rubinov  [3;  (k.6)]    for  stationary 
flows.  Let  Mc  =  -[v+1  +  v/(t+1)2+4].  Then  (3. 12)  gives  m'(0,0)  <  M  when 
K(0,0)  =  K' (0,0)  >  0  and  m' (0,0)  >  Mc  when  K(0,0)  =  K« (0,0)  <  0. 

(D)  Suppose  that  K(0,0)  <  00 .  Then,  from  (3-5),  we  have 

lim  (cos  n«)  g(s,0)  =  -2p  (3.1*0 

is  finite  and  (3  >  0.  Hence, 

cos  u'(s,0)  ■  -  1  -  2pfl  +  o(s)  ,  (3.15) 

where  o(f)  means  that 

lim  2M  =  0  . 

s*0 


Hence, 


Since 


2 
sin  |i'(s,0)  =  2ps  +  o(s) 


lim  k  sin  |i'(s,0)  =  lim  (f       -[i1      )  sin  u'(s,0)  =  lim  (cos  n)   (s,0), 

I  eZtB  ,  S  ,S 


k(s,o)  *  +  l±  +  o(^r)  , 


and  then 


a(s,0)  =  I  +  ./^  +  oGTs) 
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\2(s,0)  =  sin  a  =  sin(|  +  y20s  +  o(7s))  =  cosCyspi  +  o^))  =  1  +  o(s), 
where  0(f)  means  that  -^-     is  bounded  as   s^o„ 

^(s^O)  =  cos  a  =  sin(+>/23B  +  oG/ir))  =  7  J2&E  +  oQt)    . 


Hence, 


a   (s,0)  =  s  +  0(s2)    . 


(3.16) 


^(8,0)  =  T|^s3/2  +  o(s3/2}  =  TMs3/2  +  o(s3/2}   g  ft 


3V  K     T  "v°   '  -  *  v  — 

Therefore  the  shock  is  approximately  of  the  form 

/  9   2x1/3 
a2  =  (Bfiai  )  '     •  (3.18) 


Let  P,  =  /3*  .   Thpn  a   =  J.  A   o  3/2 


1  ^v  9  *  Then  ai  =  i  pi  V   and  from  (3°5^ 


(^  +  7  -  1)K(0,0)  =  (2m.2-3-7)K'  -=§-■»,  _  ^  +  ^^,2      ^  j 

7v   'S    v'2      *  (2+(r-l)m'2    jUn 

(3-19) 

(3.17)  and  (3.18)  extend  Lin  and  Rubinov's  result  [3; (5.5)  and  (5=6)]  to 
a  general  stationary  and  pseudo-stationary  flow,  and  their  equation  (5.6)  should  be 


(a%  -  (7  +  3))  i 


2 


/  2   ,    >j  1   Qk2  (        27MT  -  ^-1 

1   (V         J  2       (    2[l+(7-l)^] 


(3.20) 


H^  =  m',  since  in  their  case  rry—  (0,0)  =  1  . 


^'     Normal  Shock  Attached  to  a  Straight  Wall 

In  this  section  we  will  assume  that  the  flow  ahead  of  the  shock  is  uniform 
and  the  continuity  and  differentiability  requirements  on  functions  involved  below 
are  assumed  satisfied. 
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Then 


and 


where 


and 


ii(s>T)   =  J2,T(S'T)  =  °  C^-1) 

cr   _(s,0)  =  km'  sin  ji1  (2^.2) 
>b 

ji,T(s'0)-Aioej^  ».3) 

61(s,0)  *  e±  k  ,  (U.U) 

2   2    2 

e£(s,0)  =  -jj-  m  sin  n  ,  (k.5) 

e2(s,0)  -=§£&  -i~+*i2J£)  .                          (^.6) 

m    <3 

Since  the  shock  is  normal  to  a  straight  vail,  by  taking  s  =  0  to  be  the 


wall,   we  have 


al 


;(0,0)  =  ■£,        J2(0,T)  =    J2(0,T)  =  0,    or  ff.  (^7) 

For  simplicity,    let     p    (0,T)   =     J»2(0,T)   =    TT  i     Then 

u.(0,0)  =7To 
Suppose  that  (\i-lf)  Ti(0,0)   =  0,    for  i  =  0,    1,    . «.,   h-1    .     Then 

(sin  n)  t1(0,0)   =   (cos  n    »   u  T)  Ti_i  =  0  (k.8) 

and 

^Th(0,0)  =  -a^h  -  -(SiHJt  J2;S)^h.i  =  0 

by  (2.5*0  •     Therefore,    for  n  =  0,   1,    2,    3,    ■  ••, 

(n-7T)^Tn(o,o)  =  ^ja  =  A22,Tn  =  °  '  (1k9) 

From  (k.k), 

5,(0,0)   =  0      . 

Let 

51  Ti(0,0)   =  0,        for  i  =  0,   1,   2,    ...,   h-1    . 

Then 
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h      h 


0  =    j2jTh+l(0,0)  -   (A2J  8.))Th  .     Zq   (1)A2.)Tl  5.^ 

h 
=  ^(i)A21^i  61^h.i  =  A21  5^Th     . 

Since 

A21(0,0)   ^  0   ,  &^Th(0,0)   =  0    . 

Hence 

5-^Tn(0,0)   =  0        forn  =  0,    1,    2,    .„  (4.10) 

In  the  following,   we  will  show,   by  induction,   that 

5l,s2iTn    (°>0)   =  52,s2i+lTn  -  0  (4.11) 

for  i,n  =  0,   1,    2,    ...    .      Suppose  that   for  i  =  0,   1,    ...,i,    n  =  0,    1,   2,  . .., 

,   hx^nlOrO)  =  0     ,  {h,12) 

5l,s21Tn    (°'0)   =  °     '  (^.13) 

52^s2-lTn(°^°)   =  °  >  (4.l4) 

m,s2i-lTn(0,0)   =  0  ,  (4.15) 

c^s2i-lTn(0,0)   =  0  ,  (4.16) 

(n-77)^s2iTn    (0,0)      =  0  ,  (4.17) 

(sin  n)^s2iTn(0.0)   -  0     ,  (4.l8) 

k  R2i-1    (0,0)      =  0  ,  (4.19) 


and 


where  f     imj    (0,0)  =  0  if  i  or  j   is  negative.     To  complete  the  induction,   we  need  the 
following  lemma. 
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Lemma.      If  f 


j  =  0,   1,    ...,  JL  ,   then 
,-1 


,s2i+iTj   (0,0)  =  g^s2i+lTj  =0     for  I  =  0,  1,   ...,  n  aad 


(f-    )^s2i+lTJ    (0,0)   =    (fg)^21+1  Tj   =  0      . 


The  proof  follows   from 

,-1n  ,_   _x  ,1  ■  21     2i 


(f_    'i***^0'0'  *   "^  ',*>,•**  «   "  q^   (»  KAj)  rf  f  s2i+1.q]  . 


(fg),s2i+irJ  (o,o)  =  fr  a  g  +  fg^s)^2iTj  =  o  . 

By  the  lemma, 


^S2i-Lpn(0,0)   -    (mc)^2i.lTn  .  0  (4g2o) 


t        \       1   /  ,        2i2i 

<T7s2i+l    (0^0)   -  ~(kv  sin  n)^gl  =  —    £  (j    )for)     j(«tn  n)^^   =  0      (4,2lj 


(cos  n)^2i+lTn(0,0)  =  -(sin  W.)     2Lp  =  0        .  (fc.fe) 


V     P-v  .  O  /  p" 

2/  (27  O-  +I-7)        COS^    (i 


i^trrx  .^— i-  c*-23) 


mf82i+l    (0,0)  -  0  , 


^^s^s21    (0'0)   =  2m,s21+1  =  0      *  (^.2*0 

y^s2i+l    (0,0)   .   Il  +  Q-tcr2-!)]  ^2i+1  =  0    . 

2 
■tf    s2i+l    (0,0)   =    (7+I)    [—2= g]      Pi+1  =   0    . 

2  ^1 

=jS2i+l    (0,0)   .  —  (^))S2i+l  .  0   . 
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Therefore, 

■,.21+1    (0/0)  =  0.  (^>25) 


Let 


M^s2jl+lTj(0,0)   =   c^2i+lTr  0,      fbr  J  »  0,   1,    ...,    h-1. 
Then  by   (2.51)   and   (2.52) 

»  fl2*+l*h   (0,0)  .   -[(i  +  2^i  m)J^T  +  f]^^ 

2^+1       x 
=   -EC      q   H(5  +  <§1  m)^q  ]1^Ts2Jl+l-q]^1h-l 

2£+l 


and 


7-1       f  7-1        2^+1  \ 

!,S2  +lTh   (°>0)   ■  -2~(c  Ji,T),s2/+lTh-l  =     2"  S(      1   Hc,sq  01>Ts2i+l-q]  ^Th-1 


7-1 


i  2^+1 
^s2/+lTh-l  =     2"  ^< 

2£+l 


t~E(      q    )[c,Sci5l,Ts2Z-q^Th-l  =   0      ' 
By   (4.24)   and   (4.25),      for  n  =  0,    1,    2,    ..., 

^s2|+lTn(0,0)   =   c^+Lpn  =  v^+Lpn  =  0      ,  (4.26) 


and 


since 


( g  1     g     ):B21+Lpn   (0,0)  -  0  ,  (4.27) 

l-m     cos     |i   ' 

(l-m     cos     |i)   s2i+ijn   (0,0)  =  0     . 

2  21+1  2 

A12,s2i+1    (0*0)   =    (      \   COS2H),s2i+l  =  S(      q   )(-         m  )         (cos  n)      2i+l-q 

l-m     cos  (i   *  l-m  cos  jj.   '  ' 

=  0   . 

2#+l  ^ 
0  =  5^^(0,0)   -    (a^o.)^^!   =  E(      q   )A1AtSq  6^s2JUl-q  =  A^^    . 
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Hence, 

B2,g2(Ul(0,0)   .  0.  {Ks6) 

Let 

52,s2|+lTj    (°'°)   =  \s2^+V   =  °>      for  J   *  0,    1,    ...,    h-1    „ 
Then 

=2,s21+itJ  (0,0)  .    (iLSSSJLy^^  .  Vofc  ,,(304*)^.,] 


\  T(s,T)  =  cos 


^2,s      •  (U,29) 


Hence, 


^.Zt+Lph   (0,0)  =    (cos  ^^2fS)ts2UlT^l  =   [««  u(62+c2n,s2i+lTn-l 

2X+1 
=  Z(      q   )t(cos  M);Sq(52+c2)^2ifl.q]^i1.i  =  0  , 

2H+I 

=2,s2l+iTh(o,o)  =  E(    ,  )[X)Sg(H2^))s21+1.q])Th  .  0  . 

2 

A12,s2i+lTli    C°*°)   =   [m     co°  B         ■]      2i+lmh  =  0   , 

\(l-m  cos  |i)    ' 

0  =  5^s2iTh+l    (0,0)   =   (fhT)fSd+h*  «    (VP,h2^ 

2l+l 
=  E(     q   )(Alj,SCL  5j,s2i+l-q),Th  =   TA12  S2fB2l+l\Th  =  A12  62^2£+lTh   ^ 

rherefore,    for  i  =  0,   1,    2,    ..„,/      and  n  =  0,   1,  2,    ..., 

52^s2i+lTn    (0,0)   =  ^s2i+lTn  =  A12,s2i+lTn  =  C2,s2i+1T*  =  °   '    ^'3°^> 


and 


Now 


and 


iince 


re  have 


(sin  n)      2i+l    (0,0)   =  2(sin  n  cos  ^      )      2i  =  0   , 


s    ,  s' 


>°   > 
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.    2 
/Sin     |iv  ,        , 

( f)    21+1  (0,0)  -  0 


"by  the  lemma,   and  then 
b 


Hence, 


Since 


Hence, 


Therefore 


Let 


Then 


2,s2i+l    (°>°)   =  J5X&   "  ~  +  ^^  >,H2i+l  ■  0   •  (^.31) 


2J?+1 
0  =    (*V,s2J?+l    (0,0)   =  Z(     q    )k^q  b^i+l-q  =  k^+1    .  (*•») 

^s2i+2    (0,0)   =    (62+C2-k)^21+i  =  0    .  (U.33) 

2i+l 
(sin  (O^ate    (0,0)   .  Z(    q     )(cos  n)^q  ^s2^+2_q  =  cos  nn^j^  =  °< 

2/+1 

A11,S2^2    (0,0)   =  Z(     q    )(  X        ■■)   sq   (sin  »i)     2i+2       =  0   . 

1-m  cos  u   ;  > 


82;S2i+l  T(0,0)   =    (j2,s-C2)^+lT  -  J2^s2i+2T 


2JU-2 


=    (A2jV,s2i+2   .  Z(     q     )(A21^q51^2^+2_(1  +  A2^sq5^s2i+2.q) 
=  A21  5l,s2^+2    - 

5l,s2£+2    (0,0)    =   0      .  (km3k) 

6l,s2^+2TJ    (°'°)   =  ^,s2-^+2TJ   =  0,    for  J   =  0,    1,...,    h-1. 

2J?+1 
(sin  n)^2je+2rrj    (0,0)    =  Z(     q    )[(cos   »i)      q  ^g2J?+2-q]    T. 


=    (cos  M**j82A+2),Td   =  0   I 

^,S2?+2TJ  (0,0)  -  [   Sl2%   ]  s2g+2Tj  -  0  , 

\(l-m  cos  (aj 
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and 


Therefore 


Thus 


2.G+2 
-(sin  p  --p)^^^^  .  Z(     q   )[(A2^sq  5.^+2_q) 

-(sin,)^q(%)^e£+2.q]^1  =  0   . 

(sin  u)^2^+2Th    (0,°)   -  Aii,s4+2Th's  °      ■ 

0  =  52^2|+lTh+l    (0,0)  -    (j^g    -  G2)jS2J+lTh+l  =    (A2.5.)^2^+2  h 

2)1+2 
=  Z(      q   )(A21^q  5^s2  +2-q  +  A^q  5^ s2i+2-q),Th 

=  (A21  5l,s2i+2)^Th  =  A21  S^s2j+2Th   . 

5l,S2£+2Th  (0,0)  =  0  . 
By  (4.43)  and  (k.kh),   we  have  for  i  =  0,  1,  ...,//,  n  =  0,  1,  2,  , 

^,S2l+2Tn  (0>°)  "  \fS2jU2Tn   =  A11^s2i+2Tn  =  (s^  n),s2^+2Tn  =  0  .    (4.35) 

By  (4-35),  (4.32),  (4.30),  (4.26),  (4.10),  (4.9)  and  (4.8)  we  complete  the 

induction.  Therefore,  (4.11 )  and  (4.19)  are  true. 

t>    i    Bj    t_  (     r,\        4m2sin  (i  cos  u      ,     _      H     _   _ 
Remark .   Since  b  (s,0)  =  *-= —  ,  we  have  for  n  =  0,  1,  2,  ..., 

4   2 
bl  s2n+l  (0,°)  ■  ^Ti~(m  sin  »*  cos  ^)  s2n+1  *  °  > 

by  (4.22),  (4.24)  and  (4.l8).  Relation  (4.31)  says 
b2,s2n+l  (0,°)  =  0   • 


Hence, 


Therefore,   we  are  unable  to  determine  k     2n   (0,0)      from 

> s 

51?s2n(0,0)  =  (^\)9B2n   =  0  ,  and  6^s2n+l  (°,°)  =  (kV,S2n+l  =  °  > 
since  in  both  of  these  equations,  the  coefficients  of  k  2i  are  zero. 
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i.  Approximation  Forms 

In  this  section  we  assume  that  the  flow  ahead  of  the  shock  is  uniform  and  the 
lock  T  j=  0  is  normal  to  the  solid  boundary  s=0.  Thomas  [9]  proved  that  if  the  shock 
ias  finite  curvature  at  the  point  of  attachment  (0,0),  then  the  flow  is  not  stable, 
lere  we  assume  that  k(0,0)  =  00,  J  m(0>0)  <  °°>  and 

liTnf>TJ  <S'°)  =  f,siTJ(°'0)     '  (5.1) 

for  i,j  =  0  or  1,  where  f  4  j(0,0)  may  be  infinite.  Since  the  shock  is  normal, 

sin  n(0,0)  =  0  ,      a(0,0)  =  -|-  .  (5.2) 

?rom  (^.3) 

J2jT(s,0)  :=  A.,  e.  k  .  (5.3) 

k(s,0)  =  111  (l-*2  cos2  n)  K(s;0) ,    (5A) 

/  2  .  v       .2     .    /-    1    sin  [x 
(m  -ljcos  (a.  sm  |j.  +  sin  n(l ~  +  5— £- 

m    o 


ience,  „ 

(7+l)mTC(s,0)    ^  ^5<5) 

2(2m2+l) 


sin  u  k(s,0) ~ 


Since  |>;  p.  315] 

C^O—  COS    LL  _  .     2 

/       ~\  .  r    1  /       -,        2     \     Sin    Ui  G    COS    LI  /_    Cs 

^B(B,0)  =  -k[ +  (r-l-^)  -jaf]  -  — r-^  ,  (5.6) 

m   (7+1) v' 

sin  u   k(s,0)  CT P K(s,0)    .  (5-7) 

2 (2m +l)v 

I       (a,0)  «  K(s,0)  =  K  sa  ,     KQ  ^  0,   0  <  a  <  1=*0  <  5  <  1,  and  u(0,0)=X 

(5.8) 


,et 


'hen  2/      _  %., 

m   (7+1)  V  -+1 

cos  u(s,0)  ^  -1+ = ■ K„s 

2(2m  +l)(5+l)v 

2  /      ,  N     . ,  „   ^  -  Of  1 

im 


_!(7+l)    T'  K    \2 

sin  u(s,0)  ^ g 2  V  ■  2  (5.9) 

(2m  +l)(a+l)vj 
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and 


K0^° 


(5.10) 


Hence, 


fm2(7+l)KAv(at+l)Y/2       a-l/2 
k(s,0)^  °p 

k (2m +1) 


(5.H) 


From  (5.8), 


>("t: 


a 


)   =    fo(s,0)  -Kn  suT 


(5.12) 


Now 


W 


(s,0)  =  C0  +  c  k 


2  2 


2      /        (a+l)Knv 
1-m  ^        '0 


\l/2 


a-l/2 


'—      m 


(7+1) (2m +l)v' 


(5.13) 


Hence, 


J 


i(B,o).TT+2iksQ 

!  m 


Kov 


1/2 


(a+1 )  (7+1 )  ( 2m  +1 )  ccr  ^  v ' 


sa+1/2     (5.1U) 


and 


]2(s,T)   . 


m 


(a+l)  (7+1)  (2m2+l)c1cr  =^v' 


1/2        _ 


a+l/2 


K0ST    . 


x(s,T)   =  a1(s,0)   +  cos  \  dT    ~  a1(s,0)   -T 


(5.15) 
(5.16) 


a2(s,T)   =  a2(s,0)   + 


Jo    sinl 


,  dT  ^a2(s,0)   +    (7T-  )o(s,0))T„         (5.17) 
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fow 


3'     (s,0)-lT=   0(8.)       . 


'here fore  •   2/     ,  x  \l         5+1 

7f       y  ( m  (7+1)  v     v 

a(s,0)   --*■...  h  -n'  .  0(s)  -J 12— 

^  /  (2m  +l)(a+l)v« 


fence, 


m2(7+l)K0v        >3        S& 

^  i    ra-    2     ,r-     ~ /         s  '  (5'18) 

(fern  +l)(a+l)v» 


\2(s,0)  =  sin  a  ~,1 

fm2(7+l)Kv         -)2         %L 

\,(s,0)  =  cos  a^l g >■  s  ' 

I    (2m  +l)(5+l  )v' 

md 

a2(s,0)^s  (5.19) 

2        ^.m2(7+l)K0v  >)    I       ^1 

a   (s,0)/^^-     V » 2 /  s  2        .  (5.20) 

3+a      £  (2m  +l)(a+l)v' 

'hus,   the  shock  is   approximately  of  the  form 

a     fM2(7+l)K0v        ^1        *|3  3|5 

a  =  - —  / * \        a     =  pa      ,  say.         (5-21) 

1   3+5  I   (2m+l)(a+l)v'J    2       2 

It  is  easy  to  see  that,  if  a  =  0,  (5o2l)  is  the  same  as  (3.18)  obtained 

)efore.  Now  (5.16)  and  (5.17)  become 

a+3 
ax(s,T)  =  ps  2   -  T  (5.22) 

1   - 

.(.,,) . .    s!k£l  f—^L — ?'  ."r ,  .      (5.83, 

*         t   m     )   (5+1)  (7+1)  (2m  +l)v'j 
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